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Neutron diffraction has been used to investigate the magnetic correlations in single crystals of the 
spin ice materials Ho2Ti207 and Dy2Ti207 in an external magnetic field applied along either the 
[001] or [110] crystallographic directions. With the field applied along [001] a long range ordered 
groundstate is selected from the spin ice manifold. With the field applied along [110] the experiments 
show that the spin system is separated into parallel (a) and perpendicular (8) chains with respect 
to the field. This leads to partial ordering and the appearance of quasi-one-dimensional magnetic 
structures. In both field orientations this frustrated spin system is defined by the appearance of 
metastable states, magnetization plateaux and unusually slow, field regulated dynamics. 

PACS numbers: 75.25.+Z, 75.40.Cx, 75.50.-y, 75.60.Ej 



I. INTRODUCTION 

The spin ices 1 ^, such as H02T12O7 1 ' 3 , Dy 2 Ti 2 0^ and 
Bx>2Sn207 5 are frustrated Ising magnets. The rare earth 
atoms are located on a pyrochlore lattice, a face cen- 
tered cubic (fee) lattice with tetrahedral basis, creat- 
ing a network of corner- linked tetrahedra (see Fig. 
In these systems the crystal field creates a very strong 
anisotropy estimated at « 300 K 6 which constrains the 
spins to the local trigonal axes (members of the (111) 
direction set). The spins are therefore Ising-like but non- 
collinear: they point either into or out of the elementary 
tetrahedra of the pyrochlore structure 1 ^. The dominant 
interaction is the dipolar interaction which favors a lo- 
cal arrangement of spins on a tetrahedron in which two 
spins point in and two spins point ouli This is anal- 
ogous to the Bernal-Fowler "ice rule" for proton loca- 
tion in ice and means that the spin ice has the same 
configurational entropy as ice 2,4 . The anisotropy is so 
strong that the fully saturated ferromagnetic state can 
be reached only in principal, using magnetic fields well 
outside the range of conventional measurements 8 * 9 -. Most 
of the interesting physics happens in much lower fields, 
where the Zeeman energy, which competes with relatively 
weak exchange and dipolar interactions, is dominated 
by the (111) anisotropy. For instance, in polycrystalline 



FIG. 1: From top: groundstate of a single tetrahedron with 
field (shown by arrow) applied along [001]; groundstate with 
field applied along [110], the field separates the spins into a 
(parallel to field) and 8 (perpendicular to field) chains; long 
range ordered structure with 8 chains all parallel (called Q = 
0); 8 chains antiparallel (called Q = X) (both reproduced 
from Ref.Q). 



Dy2Ti207, about half of the zero-point entropy missing 
in zero field, (l/2)i?ln(3/2), is recovered already in a field 
of 0.5 TA 

In general the spin ices have strongly anisotropic bulk 
properties&SiiSi. This is unusual for cubic systems but can 
be understood as another manifestation of the strongly 
anisotropic nature of the individual rare earth moments. 
The application of magnetic fields at low temperature 
can be used to systematically investigate the removal 
of the degeneracy of the spin ice state. As a result of 
the anisotropy, different degeneracy removal schemes are 
expected depending on the direction in which the field 
is appliecUi. Several scenarios have been investigated 
by experimental 4 i 8 i 9 i 10 i 12 i 13 i 14 i 15 i 16 i 17 i 18 i 19 i 20 and theoret- 
ica l 11 ! 21 ! 22 ' 23 ! 24 ' 25 ! 26 ' 27 methods but there have so far been 
only a few neutron scattering investigationsiSS* 2 ^ .. 

The basic behavior of the system can be understood 
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FIG. 2: The phase diagram of the near neighbor spin ice 
model with field applied along [001] 11 . With the field applied 
in this direction the model shows an unusual liquid-gas type 
behavior. A line of first order phase transitions from the gas- 
like region (G) to the liquid-like region (L) is terminated by 
a critical point where the transition becomes continuous. 



in terms of the near neighbor spin ice model, which con- 
siders near neighbor ferromagnetic coupling between the 
spinsii. Application of the field along the [001] direc- 
tion is expected to completely remove the degeneracy 
that exists in zero fieldii. All the spins on a tetrahe- 
dron have the same relation to the applied field, making 
an angle of 54° to it (see Fig. 0) . Fields up to 2 T along 
[001] are insufficient to overcome the anisotropy and fully 
align the spins. However, the field does stabilize a ground 
state in which the same member of the set of six "two- 
spins-in-two-spins-out" local groundstates is selected on 
every tetrahedron. Experimentally 8 ^ this is manifested, 
at 1.6 K, by a simple magnetization curve and a satura- 
tion magnetization of 5.78 /zb atom -1 . The saturation 
magnetization is reduced from 10.0 /ib atom -1 , as ex- 
pected for the non-collinear structure. 

The near neighbor spin ice model predicts a particu- 
larly interesting phase transition when the field is applied 
along [001)ii. In zero field the model system shows no 
ordering transition and the heat capacity consists of a 
single Schottky anomaly due to isolated spin flips. In 
small fields a very sharp peak occurs below the Schot- 
tky anomaly, indicative of a first order phase transi- 
tion. The field-temperature phase diagram shown in 
Fig. [21 has a line of symmetry-sustaining first order phase 
transitions, terminated by a critical end point. This is 
analogous to the density variation seen in a pressure- 
temperature phase diagram of a liquid-gas system. The 
"liquid" and "gas" phases of the magnetic system are dis- 
tinguished thermodynamically by their degree of magne- 
tization, that plays a role analogous to the fluid density. 
A triple point occurs at H = T = 0. One aim of the work 
described in this paper is to search for the critical point 
in Ho2Ti20y. This is described in section ITlI Bl 

When the field is applied along [110] two spins per 
tetrahedron have a component in the field direction, 
while two are precisely perpendicular to the field (see Fig. 
|IJ. Of the two which have a component along the field, 
one is pinned pointing "in" and one pointing "out" of 
the tetrahedron. The perpendicular spins are not pinned 
by the field. However, in the applied field, the ground 
state of the tetrahedron can still be described by the "ice 
rule". Therefore the unpinned spins also form an in/out 
pair (see Fig. ^! 1-2 . The two groups form in-out spin 
chains, running either parallel or perpendicular to the 
field. Each chain has a magnetic moment in the sense 
of the chain, but in the near neighbor model there is no 
coupling between the perpendicular chains. We will use 
the terminology of Hiro: 18 , who named the chains parallel 
(perpendicular) to the field a (/?) chains. 



In the neutron scattering study of Harris et alX a single 
crystal of Ho2Ti2 07 was aligned with the field applied 
along [110]. In zero field Ho2Ti 2 07 has a disordered ice- 
like state with zero point entropy^i^. In Ref. Q it was 
shown that the application of the field along [110] restores 
order, with two coexisting magnetic structures, denoted 
Q = and Q = X (see Fig. 0. The structures are 
named because Q = gives rise only to magnetic Bragg 
scattering at zone centers, or Q = positions, whereas 
Q = X gives rise to Bragg scattering at both the zone 
centers and the X points on the zone boundaries (see 
Fig. |3J). The Q = structure was suggested to have 
all the chains, both a and (3, aligned parallel with their 
neighbors. The moment of the structure is at 45° to 
the applied field. The Q = X structure was suggested to 
have alternate (3 chains antiparallel, with the net moment 
in the direction of the applied field. It was found that 
the Q = structure was formed at 0.3 K when a field 
was applied and the Q = X structure appeared as a high 
temperature (« 1 K) modification. The Ho 3+ moment 
is severely reduced in the Q = X phasei. One of the 
results of the current work is to show that the Q = and 
Q = X structures are in fact associated with different 
sets of spins: the a and [3 chains respectively and that 
the formation of Q = X correlations is controlled by 
the very slow dynamics in the T < 0.5 K regime (see 
sections IV Al and IV B|) . 

In addition to making the first determination of the 
residual entropy in a spin ice^, Ramirez and co-workers 
measured the heat capacity of a Dy2Ti20y powder in 
fields up to 6 T. Three small, sharp, field independent 
peaks were observed in finite fields. These peaks at 
T p = 0.34, 0.47 and 1.12 K were attributed to the or- 
dering of spins not pinned by the field. More recently, 
single crystal samples of Dy2Ti2 0y have been studied by 
bulk methods in fields applied along [110] 18 . The exper- 
imental results of Ref. suggest the formation of the 
[3 chains, but do not infer any correlation between the 
chains. It was concluded that the (3 chains behave as one 
dimensional Ising ferromagnets. Most recently, simula- 
tions of the dipolar spin ice model (that goes beyond the 
near neighbor model by considering further range dipo- 
lar couplings) have predicted that in strong fields there 
are two transitions 24 - 25 . At high temperature the field 
induces order in the a chains and at a lower temperature 
the (3 chains undergo a phase transition into the Q = X 
phase. 

Dy2Ti207 has been studied less by neutron scattering 
due to the significant absorption problem. However, it 
has been shown that at high temperatures (T > 4 K) 
Dy2Ti2 07 is paramagnetic with fluctuations dominated 
by single ion effects^. A preliminary report of some of 
the results detailed below has also been madeS 8 .. 

It is now generally accepted that the spin ice materials 
Ho2Ti207 and Dy2Ti207 only differ in detail; a compari- 
son of the two should identify generic behavior and reveal 
the importance of minor terms in the spin Hamiltonian. 
In this work Ho2Ti207 and Dy2Ti207 have been stud- 
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ied by single crystal neutron scattering with an external 
magnetic field applied along the [001] or [110] axis. These 
measurements are the first detailed microscopic investi- 
gation of the field dependent properties of Dy2Ti2C>7 at 
low temperature, and an extension of the investigation 
already performed on Ek^T^O-ii 

The rest of the paper is organized as follows. Following 
an experimental section lpT|l . results are given for applica- 
tion of the field along [001] in section lTTll and discussed in 
section Hvl Then the results for the [110] orientation are 
given and discussed in sections lyland lvTl As yet there 
is no detailed theory with which to analyze these results, 
so we restrict ourselves to largely qualitative discussion 
and conclusions fsection IVII|> . 

II. EXPERIMENTAL 

The experiments on Ho2Ti2C>7 and Dy2Ti2C>7 were 
conducted on different instruments at two neutron 
sources: Ho2Ti 2 07 on D10 at the ILL reactor source, and 
Dy 2 Ti 2 7 on PRISMA at the ISIS spallation source. 

All the measurements reported in this article used a 
single crystal aligned with either the [001] or [110] axis 
parallel to the applied field (vertical) . This allowed us to 
access the scattering plane with wavevectors of the type 
h, k, or h, h, I respectively (see Fig.|21for an illustration 
of these planes) . 



A. Dy 2 Ti 2 7 on PRISMA 

Two experiments were carried out on the PRISMA 
spectrometer at ISIS under the same conditions. To mini- 
mize neutron absorption, the isotopically enriched single 
crystal of 162 Dy2Ti2C>7, described in detail in Ref. l3ll 
was used. The crystal was varnished into a large copper 
support to ensure good thermal contact and mounted on 
a dilution refrigerator insert in a 7 T cryomagnet. 

The 16 detectors of PRISMA make radial trajectories 
across the scattering planed. In order to map the scat- 
tering plane, the crystal is rotated about the vertical axis 
and scattering from adjacent sectors of the scattering 
plane is collected in successive settings. To cover a small 
area of the scattering plane two different types of scan 
can be used. Either a single setting of the crystal provid- 
ing relatively coarse coverage of reciprocal space, or for 
better resolution a rocking scan. In this case the 16 de- 
tector trajectories of successive points in the rocking scan 
are interleaved, considerably improving the definition of 
features in reciprocal space. The scattering planes and 
detector trajectories used in this study are illustrated in 
Fig. El 

Where non-resolution limited features were observed 
they were used to extract correlation lengths. The reso- 
lution function as a function of time of flight (parallel to 
the detector trajectory) is a bilateral exponential func- 
tion describing the pulse shape from the moderator™. 



Perpendicular to the detector trajectory the resolution 
function is Gaussian. The resolution functions were pa- 
rameterized and interpolated by fitting to Bragg peaks 
in the high field data and convolved with a Lorentizian 
to fit the experimental data. 



B. Ho 2 Ti 2 7 on D10 

Ho2Ti207 was measured on the single crystal diffrac- 
tometer D10 at the ILL. Two single crystals grown from 
a flux of lead fluoride were used^i, one for each field 
direction. For the [001] orientation the crystal was 
« 3 x 3 x 3 mm and for the [110] orientation the crys- 
tal was « 2 x 2 x 2 mm. They were of identical quality 
if the widths of Bragg peaks were compared. To sim- 
plify corrections for effects such as demagnetization and 
extinction, the crystals were sphericalized^. 

Appropriately aligned crystals were fixed to oxygen 
free copper pins using stycast resin and attached to a 
dilution refrigerator insert in a 2.5 T vertical field cry- 
omagnet, mounted on D10 using a tilt stage. D10 was 
configured with a graphite monochromator and a small 
area detector. The wavelength was 2.56 A and a py- 
rolitic graphite filter in the incident beam was used to 
suppress any A/2 contamination, oj-scans (in which the 
crystal is rotated about the vertical axis) were measured 
at all positions. 

The integrated intensities of nuclear and magnetic 
Bragg peaks were extracted in the usual manner 36 . How- 
ever, the diffuse scattering was often too broad for the u>- 
scan performed. This resulted in the appearance of a flat 
intensity profile, albeit elevated from the expected back- 
ground, whose intensity and width were subsequently 
strongly field or temperature dependent. It was therefore 
attributed to a magnetic scattering feature broader than 
the width of the scan. It proved impractical to extract 
widths from these features so the peak intensity was ob- 
tained to give a measure of xq, the wavevector dependent 
susceptibility, throughout the experiment. The ampli- 
tude was extracted by fitting to a functional form appro- 
priate to the peak shape: straight lines, single Gaussians, 
or triangular functions. 

The Bragg scattering data was treated by refining 
model magnetic structures using the CCSL 37 program 
MAGLSQ. During the data analysis it became appar- 
ent that there was a significant extinction problem. In 
a subsequent experiment the larger crystal was mea- 
sured at room temperature using four circle geometry and 
two wavelengths (1.26 and 2.36 A, copper and graphite 
monochromators respectively). This data was used for 
the refinement of extinction parameters using the Becker 
and Coppens model as implemented in the CCSL pro- 
gram SFLSQ 37,3 —. Thermal parameters derived from a 
Rietveld refinement against powder diffraction data were 
used to avoid correlation between the extinction and ther- 
mal parameters^. The powder data was recorded at 1.8 
K on the POLARIS medium resolution diffractometer at 
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FIG. 3: The reciprocal space of Dy 2 Ti 2 7 in the h,k,0 (left) and ft, h,l (right) planes as studied on PRISMA. The Q = 
and Q = X positions are indicated by o and x respectively. In both diagrams the area of scattering plane that lies within 
the first aluminium powder line and was included in a mapping is indicated by the curving shells. With the field along [001] 
(ft, fc, scattering plane), single setting scans were measured at intermediate field points. In this case 16 detectors cover the 
region between the two radial detector trajectories enclosing the ft, 0, axis. With the field along [110] (ft, ft, I scattering plane), 
rocking scans were measured. In this case 80 detector trajectories cover the region between the two radial detector trajectories 
enclosing the 0, 0, 1 axis. This gives vastly improved reciprocal space coverage compared to a single setting scan such as that 
indicated about the ft, ft, axis, where only 16 detector trajectories cover the same area of reciprocal space. 



ISIS. The difference in temperature of the single crys- 
tal and powder data sets was assumed to be absorbed 
by the extinction parameters as the final fit was good 
(i?3 = 4%2 9 .). The extinction parameters allow the re- 
production of nuclear intensities so that the magnetic 
structure can be fitted by adjusting only the size of the 
magnetic moment. 

The quality of the fit of the magnetic structures de- 
creases as the field increases (in Fig. it can be seen 
that as the field increases, the uncertainty of the fitted 
moment becomes much greater). This is attributed to a 
field dependent extinction effect. It is possible for nuclear 
and magnetic scattering to be affected differently by ex- 
tinction^. Since the magnetic structures in this work are 
field dependent, it is likely that the extinction is also field 
dependent. As the field and magnetic intensity increases, 
the extinction correction becomes less and less applica- 
ble. It proved impossible to overcome the problem with 
the available data. 



III. RESULTS FOR APPLICATION OF FIELD 
ALONG [001] 

A. Dy 2 Ti 2 7 

The sample was cooled to 0.05 K and the h, k, scat- 
tering plane was mapped. Then the first quadrant of a 
hysteresis loop was measured (up to 2 T), during which 
the scattering plane was mapped at 1 and 2 T. At inter- 
mediate fields, a single setting was measured giving ac- 
cess to a small area of reciprocal space surrounding the 
h, 0, axis out to the (8, 0, 0) reflection. After returning 
to zero field the crystal was rotated to access the (4, 2, 0) 
peak and its field dependence was measured. 

Fig- El compares the scattering in and 1 T from the 
h, k, plane. The diffuse scattering entirely disappears 
and is replaced by magnetic Bragg scattering at Q = 
positions consistent with the formation of the expected 
canted ferromagnetic phase. This is direct confirmation 
of the complete removal of degeneracy and field induced 
ordering. 

The formation of the ordered phase is also notable. 
Fig. [S] shows the integrated intensity of two purely mag- 
netic Bragg peaks. This magnetic intensity grows in a 
series of steps. Successive measurements produce simi- 
lar, but not identical, order parameter curves indicating 
a strong history dependence. 



FIG. 4: (Color Online) Dy 2 Ti 2 7 (B // [001]): scattering in- 
tensity in the ft, k, plane recorded on PRISMA at 0.05 K in 
zero field (top), and in 1 T (bottom). The strong diffuse scat- 
tering characteristic of the short range ordered spin ice state 
disappears entirely, giving way to magnetic Bragg scattering. 
Note that the intensity scale is cut off at 2 to reveal the pres- 
ence/absence of diffuse scattering, the maximum intensity of 
a Bragg peak is much greater. 



FIG. 5: Dy 2 Ti 2 7 (B // [001]): Integrated intensity of the 
(2, 0, 0) peak on the first field cycle and the (4, 2, 0) on the 
second cycle with field applied along [001] (the intensity of 
(2,0,0) has been divided by 2). Both peaks are purely mag- 
netic scattering. The history dependence and existence of 
metastable states during the magnetization process is clearly 
seen by the different hysteresis in the two cycles and the 
plateaux. Lines are to guide the eye only, filled symbols indi- 
cate the falling field leg. 



B. Ho 2 Ti 2 7 

In an attempt to find the anticipated liquid-gas critical 
poinljii a detailed survey of -B/T-space was performed. 
Field cycles up to 2.5 T were performed at 0.05, 0.15, 
0.3, 0.5, 0.7 and 1.2 K. Each was separated by warming 
to 1.2 K and cooling in zero field. For this study four 
reciprocal lattice positions were measured, two Q = 
positions ((2, 0, 0) and (4, 2, 0)) and two Q = X positions 
((3,0,0) and (1,0,0)). 

While cycling the field, magnetic scattering was found 
to develop only on the Q = Bragg positions. The loops 
at all temperatures are consistent with the development 
and saturation of the Q = structure (Fig. |SJ. In all 
cases i?3 < 10% 39 . The value of the ordered moment 
obtained (10.3 ± 0.4 {1b atom -1 ) agrees well with that 
expected for Ho 3+ (10.0 /zb atom -1 ). Again it is clear 
that the degeneracy is completely removed by the appli- 
cation of the field in this direction. Increasing hysteresis 
and persistence of metastable states was observed as the 
temperature was decreased below 0.5 K. At 0.05 K a rem- 
nant magnetization is frozen into the sample on return 
to zero field. 

Although no Bragg scattering was observed at the 
Q = X positions, strong diffuse scattering associated 
with the spin ice state 3 - was observed in zero field. The 
intensity of this scattering during the field cycle at 0.05 
K is compared to that of the (2,0,0) in Fig. □ The 



5 



FIG. 6: Ho 2 Ti 2 7 (B // [001]): development of Q = mag- 
netization at different temperatures, with field applied along 
[001]. As the temperature decreases increasing hysteresis and 
metastable states appear in the magnetization, although the 
final ordered state is always the same. Lines are to guide the 
eye only, filled symbols indicate the falling field leg, successive 
higher temperature loops are shifted by -1 T and + 0.5 fiB 
atom -1 . 



FIG. 7: Ho 2 Ti 2 7 (B // [001]): the intensity of the diffuse 
scattering, at the (1,0,0) position as a function of field, at 
0.05 K (top). The (2,0,0) peak intensity is shown as a com- 
parison (it has been divided by 50). As the fully ordered 
structure develops the sample passes through a partially or- 
dered regime in which both Bragg scattering and diffuse scat- 
tering are observed (between 0.3 and 0.7 T). At the end of 
the field cycle the sample is frozen in a partially ordered state 
with a remnant magnetization. As the temperature is raised 
above 0.5 K, the remnant magnetization decays and the dif- 
fuse scattering characteristic of the zero field spin ice state is 
reestablished (bottom) . Lines are to guide the eye only, filled 
symbols indicate the falling field leg. 

diffuse scattering persists in finite field and shows a step- 
like decrease commensurate with the step-like increase 
in the (2,0,0) scattering. This indicates that disordered 
regions persist in the sample in finite field. The diffuse 
scattering disappears as the Bragg scattering saturates in 
the fully ordered state. This relationship between Q = 
Bragg and Q — X diffuse scattering is true at all temper- 
atures measured, but at elevated temperatures (T > 0.5 
K) there is no hysteresis. 

After the hysteresis loop at 0.05 K was completed, the 
temperature was raised in zero field. As discussed above 
and shown in Fig. there is a remnant magnetization at 
the point where the hysteresis loop returns to zero field. 
By comparing the Bragg and diffuse scattering intensities 
during the temperature rise, it can be seen that the frozen 
long range order created by applying high fields at low 
temperatures decays and the short range ordered spin ice 
state is re-established as the system is warmed up (see 
Fig. 0). This occurs at T w 0.5 K. 

IV. DISCUSSION OF RESULTS WITH FIELD 
APPLIED ALONG [001] 

Both materials show similar behavior: they are driven 
into the Q = structure by the application of the field 
and are strongly hysteretic/history dependent at very low 
temperature. 

When measuring Ho2GaSb0 7 , Blote et al. noted that 
below T » 0.3 K it appeared that the spin system was ap- 
parently no longer in good thermal contact with the ther- 
mometes^i. Similar effects were observed by Orendac in 
more recent heat capacity measurements on Ho2Ti20 7 , 
below a temperature of T w 0.6 K the heat capacity be- 
came time dependent^. The neutron scattering data is 



comparable with this as the hysteresis develops in field 
scans at 0.5 K and below (see Fig.[6J|. The spin dynamics 
of Ho2Ti 2 07 are slow in this temperature range: no dy- 
namics have been detected by inelastic neutron scatter- 
ingH and Matsuhira et al. found that the susceptibility 
of Ho2Ti20*7 was virtually zero below 0.5 K^4. How- 
ever, Neutron Spin Echc^, and muon spin relaxations* 
measurements have observed some residual spin dynam- 
ics at these low temperatures. The same is also true of 
Dy 2 Ti 2 Q 7 9 i 47 i 48 . Below T « 0.5 K there must be a cross 
over from a regime in which the spin interactions main- 
tain the system in a spin ice state but there are sufficient 
dynamics to allow relaxation, to one which is completely 
frozen. This was clearly demonstrated in Ho2Ti20 7 by 
the decay of the remnant magnetization into the disor- 
dered spin ice state as the sample was heated above 0.5 
K (see Fig. |7J). 

The presence of magnetization plateaux in the hys- 
teresis loops at 0.05 K is particularly interesting. The 
temperature (in Ho2Ti20 7 ) and history (in Dy2Ti207) 
dependence of the hysteresis suggests that these are due 
to very slow dynamics. The system is almost completely 
frozen and must be driven into the ordered state by the 
field. An alternative way of looking at this is to ask 
why a plateau occurs at a particular value of the mag- 
netization? For example in Ho2Ti20 7 it is at w 6 /ib 
atom -1 . We propose that because the dynamics are so 
slow and the thermal fluctuations so weak, that the sys- 
tem is barely able to relax. Therefore, the plateau does 
not represent a partially ordered structure with a mag- 
netization of 6/10, or fluctuating moments contributing 
an average of 6/10 of the saturated magnetization. In- 
stead, it represents a partially ordered spin ice, in which 
6/10 of the moments have been driven into their ordered 
orientation by the field and pinned. The remaining 4/10 
require higher fields to reorient them. 

To extend this conjecture we suggest that there must 
be loops, chains and clusters of spins that can be ori- 
ented relatively easily. Meanwhile some spins are left on 
branches, or in regions not reached by a percolating clus- 
ter. In the absence of thermal fluctuations these spins be- 
come trapped as defects in the Q = state and only high 
fields can remove them. This is similar to work hardening 
in alloys. When alloys are first worked, defects can move 
to release the stress. As the working continues, remain- 
ing defects become entangled or locked together. Even- 
tually the material no longer responds to the stress and 
it fails. The neutron scattering results from Ho2Ti20 7 
support this hypothesis. The system begins in a disor- 
dered spin ice state, with diffuse scattering and no Bragg 
peaks. Applying an external field results in magnetic 
Bragg scattering from an ordered structure. Simultane- 
ously the diffuse scattering from short range correlations 
decreases. However, in the intermediate plateau region 
both Bragg and diffuse scattering are present, indicating 
that the system has regions of both long and short range 
order. 

The postulated liquid-gas critical point was not ir- 
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refutably absent. The prediction of Harris et al. can 
be cast as a first order phase transition in finite fielcUi. 
fn both materials this is observed at 0.05 K, but it is ab- 
sent above 0.15 K in Ho2Ti207. The critical point must 
lie in this temperature range, if it does exist. Using the 
effective nearest neighbor coupling strength (J e g = 1.8 
Ki£) and expected rare earth moment sizes (10.0 /j,b 
atom -1 ) the critical point is expected at T c = 0.29 K 
and B c = 0.27 T for Ho 2 Ti 2 7 . The development of the 
hysteresis in Ho 2 Ti 2 7 occurs in the right temperature 
and field regime but does not agree in detail with these 
predictions. The liquid-gas critical point is a property 
of a thoroughly equilibrated simulated system: it may 
very well be obscured in the real materials by the lack 
of dynamics in the temperature range of interest. Also, 
the behavior of the dipolar spin ice model in applied field 
is largely unknown, it is possible that the critical point 
may appear elsewhere in (B, T) space for this model. 



V. RESULTS FOR APPLICATION OF FIELD 
ALONG [110] 

A. Dy 2 Ti 2 7 

Because two experiments were performed and the path 
through (B,T) space became quite complex, we discuss 
the experiments in terms of similar field/temperature sec- 
tions. For data points between the large h, h, I maps, 
rocking curves were measured about the 0, 0, 1 axis out 
to (0,0,10). 



1. Hysteresis loops at 0.05, 0.3 and 1.9 K 

During the two experiments two hysteresis loops were 
measured at both 0.05 (one in each experiment) and 0.3 
K (both in the second experiment), and one at 1.9 K (sec- 
ond experiment). Information about the development of 
spin correlations can be extracted from the behavior of 
both Bragg and diffuse scattering. In zero field at 0.05 
and 0.3 K there is diffuse magnetic scattering typical of 
the spin ice ground stated. On applying a magnetic 
field this diffuse scattering disappears and is replaced by 
diffuse scattering of a new form centered around antifer- 
romagnetic Q = X positions such as (0, 0, 2n + 1). Con- 
comitantly, magnetic Bragg scattering appears at Q = 
positions like (0,0,2) (see Fig. [SJ - 

The new diffuse scattering is sharper than the diffuse 
scattering observed in the spin ice state and it is dis- 
tinctively distributed. To distinguish it from the Bragg 
Peaks and spin ice diffuse scattering it will forthwith be 
referred to as a Q = X feature. The Q — X features are 
relatively narrow parallel to the h, h, direction, but re- 
main quite wide along 0, 0, 1. The distinctive distribution 
of scattering in the Q — X features is clearly seen in the 
mapping of the scattering plane in 1.5 T during the first 



FIG. 8: (Color Online) Dy 2 Ti 2 7 (B // [110]): Scattering 
intensity in the h, h, I plane recorded on PRISMA at 0.05 K 
in zero field (top), and in 1.5 T (bottom). In zero field the 
diffuse scattering is characteristic of the short range ordered 
spin ice state--. In 1.5 T both Bragg (at Q = positions such 
as (0, 0, 2)) and diffuse (at Q = X positions such as (0, 0, 3)) 
magnetic scattering is observed indicating the coexistence of 
long and short range magnetic order. Note that the intensity 
scale is cut off at 2 to reveal the presence/absence of diffuse 
scattering, the maximum intensity of a Bragg peak is much 
greater. 



hysteresis loop at 0.05 K (see Fig. |8J). Even in a field of 
1.5 T the Q — X features are not resolution limited. 

At 0.3 K the disappearance of the spin ice diffuse scat- 
tering and development of the Q — X features is similar 
to that described for 0.05 K. The form of the hystere- 
sis loops differ in detail, but both Bragg scattering and 
Q = X features were observed. At 1.9 K only Bragg scat- 
tering was observed, the Q = X features did not appear. 

Structure factor calculations show that the Q = 
Bragg scattering is due to the a chains and the Q = X 
features are due to the (3 chains (see section I for def- 
inition of a and (3 chains). The shape of the Q = X 
features gives considerable information about the spin 
correlations induced in the f3 chains by the application 
of the field. The position of the features (i.e. at Q = X) 
and the fact that they are discrete in reciprocal space 
shows that there must be antiferromagnetic correlations 
between (3 chains. The anisotropic shape of the features 
shows that the order is longer ranged parallel to h, h, 
(intrachain correlation length) than parallel to 0, 0, 1 (in- 
terchain correlation length). A single correlated region 
is therefore a collection of a few chains (with an inter- 
chain corr elation length C = 14.2 ±0.7 A this is 1.13 unit 
cells (in the 1,1,0 direction) and so 4-5 chains), each ex- 
tending 770 ± 40 A (intrachain correlation length) and 
ordered antiparallel. The actual correlation lengths mea- 
sured depend on the history of the sample (these were 
obtained at 0.05 K and 1.5 T), but in the most corre- 
lated states observed they are always of this order. We 
do not know what happens out of the scattering plane: 
the Q = X features could be rods or ellipsoidal features 
intersecting the scattering plane. Rods would indicate a 
two-dimensional ordering of (3 chains in planes parallel to 
the scattering plane. Three dimensional ellipsoids would 
indicate three dimensional domains. On the balance of 
evidence it seems likely that the interchain correlations 
are three dimensional as opposed to two dimensional. 

Comparison of the development of the Bragg scattering 
during hysteresis loops measured at the same tempera- 
ture reveals that the form of the curve depends on the 
field increments, or ramp rate, used to measure the loop. 
At both 0.05 and 0.3 K, two different hysteresis loops 
have been measured and shown in Fig. [5| At each tem- 
perature one loop was recorded with coarse field incre- 
ments/rapid change in field. These loops contain steps 
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FIG. 9: Dy 2 Ti 2 7 (B // [110]): Integrated intensity of 
(0,0,2) with field applied along [110] at 0.05 K and differ- 
ent sweep rates (top) and integrated intensity of (0, 0, 2) with 
field applied along [110] at 0.05 K and 0.3 K with slow field 
sweeping (bottom, the second loop at 0.3 K was not com- 
pleted). In the first loop the field stepping regime was much 
coarser than the second (the loops at 0.05 K were recorded 
in separate experiments and the intensities from the first are 
scaled to compare with the second). At similar sweep rates 
there is a change in dynamical regime between 0.3 and 0.05 
K as shown by the presence/absence of steps at 0.05/0.3 K. 
Lines are to guide the eye only, filled symbols indicate the 
falling field leg. 



FIG. 10: Dy 2 Ti 2 7 (B // [110]): the amplitude of the Q = X 
feature at (0, 0, 3) during hysteresis loops at 0.05 K and 0.3 
K in the second experiment. The degree of correlation in the 
[3 chains is temperature dependent. At 0.3 K they become 
more strongly correlated. Lines are to guide the eye only, 
filled symbols indicate the falling field leg. 



or discontinuities in the development of the order pa- 
rameter. Subsequent investigations, using finer field in- 
crements generally reduced the number of steps in the 
loop but those that remained were more clearly resolved. 
Also, with fine field increments, there is a clear differ- 
ence between 0.05 and 0.3 K. With fine field increments 
at 0.3 K the steps disappeared and a very smooth in- 
crease in intensity was observed. At 1.9 K no hysteresis 
was observed. Again it appears that there is a change 
in dynamical regime at low temperature as shown by the 
presence/absence of steps at 0.05/0.3 K. The tempera- 
ture of the change is lower in Dy2Ti20 7 than Ho2Ti207 
(a similar effect was observed at 0.5 K in Hx>2Ti207 and 
discussed in sections IlII Bl and llVfl because the coupling 
strengths are somewhat lower in Dy 2 Ti 2 7 . As with the 
field applied along [001] the very slow dynamics of spin 
ice materials in this low temperature regime are mani- 
fested in these observations. 

The identical saturation intensity of the Q = Bragg 
peaks shows that the same long range ordered state in 
the a chains was realized in the hysteresis loops at 0.05, 
0.3 and 1.9 K. Unlike the Bragg peaks, the maximum in- 
tensity observed in the Q = X features was temperature 
dependent. Comparing hysteresis loops within the sec- 
ond experiment shows that the Q = X features become 
more intense at 0.3 K than at 0.05 K (see Fig. I10[l . At 
1.9 K the Q = X features do not form. This shows that 
for T > J off (1.1 K for Dy 2 Ti 2 7 ) the formation of the 
a chains does not require the formation of the f3 chains. 
However, for T < J e g the /3 chains form and the degree 
of correlation between them depends on the dynamics. 



2. Temperature scans in applied field 

At 1.5 T, after the first field cycle, the temperature 
was raised to w 1 K. In this temperature range, the 
Q = Bragg peaks were found to be temperature in- 
dependent. However the intensity of the Q = X features 
passes through a maximum at 0.7 K. Another tempera- 
ture scan was performed after the second hysteresis loop 
at 0.05 K. This was done in 0.4 T, to begin in the partially 
ordered plateau phase. A more pronounced maximum in 
the intensity of the Q = X features was observed at 0.6 
K, above which both Q — and Q = X scattering de- 
creased. 

In Dy2Ti20 7 1.5 T is sufficient to maintain the a chains 
in a fully polarized state, with no thermal randomization. 
The correlations in the f3 chains increase, but begin to 
decrease above a temperature of rs 0.7 K. The maximum 
temperature was not sufficiently high to leave the spin ice 
regime of Dy 2 Ti 2 7 and see complete collapse of these 
features. The second temperature scan in 0.4 T began 
in one of the magnetization step phases. Here we see 
that the rising temperature activated the development 
of further correlations in the /? chains, before thermal 
disorder began to reduce the correlations of both a and 
f3 chains. 

B. Ho 2 Ti 2 7 

The sample was cooled to 0.05 K in zero field and a 
full, 5 quadrant, hysteresis loop was measured (up to 2.5, 
down to -2.5, then back up to 2.5 T). The temperature 
was then raised to 1.2 K and then cooled to 0.05 K whilst 
maintaining a field of 2.5 T. At each field/temperature 
point at least five reflections were measured using u>- 
scans, these being (1, 1,0), (0, 0, 1), (1, 1, 1), (1, 1, 3) and 
(0,0,2). In order to determine magnetic structures, at 
certain field-temperature points an expanded collection 
of peaks was measured. Where magnetic structures were 
refined, crystal structure and extinction parameters were 
assumed to be the same as those determined for the other 
crystal used with the field applied along [001]. This as- 
sumption resulted in a relatively poor refinement of mag- 
netic structure in this experiment. We will therefore dis- 
cuss most of the results in this section in terms of inte- 
grated intensities rather than sublattice magnetization. 

1. Hysteresis loop at 0.05 K 

The entire hysteresis loop is shown for the (0, 0, 2) 
peak in Fig. ^2 The scattering at the Q = X positions 
(1, 1, 0) and (0, 0, 1) is diffuse. During the hysteresis loop 
these features never evolved into a peak-shape within the 
bounds of the w-scan. The intensity at the (0, 0, 1) posi- 
tion is also shown in Fig. 1111 

The (0, 0, 2) Bragg intensity saturates identically on 
both sides of the loop. This implies that the same long 
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FIG. 11: Ho 2 Ti 2 7 (B // [110]): Hysteresis loop at 0.05 K. 
The amplitude of the diffuse feature at (0, 0, 1) was fitted as a 
fiat line (top). The (0,0,2) peak (bottom) is directly related 
to the order in the a chains. As the a chains are ordered by 
the field, correlations are produced in the /3 chains. The /3 
chains evolve toward an increasingly correlated state but are 
only able to do so when the field produces dynamics in the a 
chains. Lines are to guide the eye only, filled symbols indicate 
the falling field leg. 



range ordered state is reached in high field, regardless 
of history and coexisting short range correlations. After 
ramping the field and driving the a spins out of the dis- 
ordered spin ice state for the first time (pristine curve) 
subsequent field sweeps are essentially all the same. The 
pristine curve lies outside the subsequent field cycles, in- 
dicating that the spin ice state is harder to magnetize 
than the subsequent state composed of a chains. 

When the Q = Bragg peaks appear, so do the Q = X 
features. Although the detailed shape of the peak is not 
accessible in a single w-scan, some further details about 
their form were established. It is possible that the fea- 
tures could be part of extended sheets of scattering. With 
the field applied along [110] the j3 chains are formed along 
[110]. A sheet of scattering might be expected to extend 
along 0,0, 1, w-scans at an (0,0, 1) position should scan 
across this sheet. Our collection of peaks allows us to 
check this possibility. Firstly, the fact that the (0, 0, 1) 
feature is flat, or outside the width of the w-scan, means 
that any such sheet is also rather broad perpendicular to 
0, 0, 1. It would be expected as a flat background con- 
tribution at (0, 0, 2) as well. The form factor of Ho 3+ 
will reduce the intensity, but both positions are quite 
close to the maximum in the form factor, so the intensity 
maximum for (0, 0, 1) in the hysteresis loop (600 counts) 
should also be clearly visible at (0,0,2). Likewise the 
presence of (1, 1, 0) as a flat feature would imply another 
sheet extending along 1,1, 1 which would be visible in 
the background of (1, 1, 1). However, the backgrounds of 
both Bragg peaks are field-independent. This rules out 
the possibility of extended sheets of scattering parallel 
to either 0, 0, 1 or h, h, and suggests that the Q = X 
features are localized, as in Dy2Ti2C>7. 

The evolution of the Q = X features throughout the 
hysteresis loop is very interesting (see Fig. lll|) . In general 
the intensity increases from the beginning to the end of 
the field cycling (the full 5 quadrants). However, it does 
so in a discontinuous, ratchet-like manner. When the in- 
tensity of the Q = peaks is saturated, there is no change 
of intensity at the Q = X positions. Only when the hys- 
teresis loop passes through a region where there are field 
induced changes in the Q = intensity do changes occur 
in the Q = X features. 

Despite the (3 chains being decoupled from the field, 
it is clear that at this temperature the field regulates 
the dynamics of both a and (3 chains. As discussed in 
section IIVI it seems certain that a spin ice is frozen as 



regards spin flips in zero field: any significant spin flip 
dynamics must therefore be produced by the action of 
the field. During the hysteresis loop it can be seen that 
not only does the field reorient the a chains, but in doing 
so propagates dynamics into the (3 chains. During these 
parts of the hysteresis loop the degree of correlation in- 
creases in the (3 chains. When the Q = structure does 
not change, neither does the Q = X, 

Two possible magnetic structures, with Q — and 
Q = X were discussed in the earlier work of Harris et 
ali. These structures were assumed to involve all four 
spins of the tetrahedral basis, which in light of these new 
results, is clearly not appropriate. We therefore used a 
magnetic structure with zero moments on the f3 chains 
and ferromagnetic a chains to fit the magnetic Bragg 
peaks (ignoring the diffuse scattering). The fitted mo- 
ment (6.5 ± 0.3 hb atom -1 ) was less than the expected 
10.0 /is atom -1 for Ho 3+ , consistent with the findings of 
Harris et alX who refined the moment on all four spins, 
assuming a fully ordered spin system. We believe the ob- 
served reduced moment is an artifact of a poor correction 
for extinction. 



2. Temperature Scans in applied field 

After the hysteresis cycle returned to 2.5 T the system 
was warmed to 1.2 K and then cooled back to 0.05 K in 
a field of 2.5 T. An enormous gain in intensity is seen 
at the Q — X positions, in particular (0,0,1), above 
0.3 K. The broad feature outside the width of the w-scan 
at low temperature, develops into a quite sharp peak. 
The correlations between the (3 chains are definitely three 
dimensional at this stage. However, these features are 
not resolution limited and do not have the form of a 
Bragg peak in the detector volume. The temperature 
dependence is shown in Fig. 1121 There is a tiny increase in 
the Q = Bragg peak intensity during the temperature 
cycle. 

Clearly the rising temperature activates the dynam- 
ics of the (3 chains and because T < J c g the restoration 
of dynamics increases correlations rather than destroy- 
ing them. During the temperature rise \Q begins to fall 
again above Trj 1 K. This is roughly consistent with the 
peak in the susceptibility seen at around 1 Ka*. Lower- 
ing the temperature again then freezes these correlations. 
Activating the dynamics of the (3 chains also allows a fur- 
ther annealing of the a chains causing the slight increase 
in Bragg intensity. 



VI. DISCUSSION OF RESULTS WITH FIELD 
APPLIED ALONG [110] 

The field-dependent properties of Dy2Ti207 and 
Ho2Ti207 with field applied along [110] show many 
generic features. In this section we discuss these and 
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FIG. 12: Ho 2 Ti 2 7 (B // [110]): temperature scan in a field 
of 2.5 T. The amplitude of the diffuse feature at (0, 0, 1) was 
fitted as a flat line (0.05-0.3 K), single gaussian (0.3-0.6 K) 
or triangular peak (0.6-1.2 K and all points on the falling 
T scan). For T < J e g increasing the temperature activates 
the dynamics of the j3 chains and also causes their level of 
correlation to increase. Once T > J e g the activated dynamics 
begin to destroy the correlations, this is evident in the down 
turn of (0,0, 1) intensity at 1.1 K. Lines are to guide the eye 
only, filled symbols indicate the falling temperature leg. 



the conclusions which can generally be drawn for spin 
ice materials with a field along [110]. 

The spin system is separated into two sub-systems, the 
a and f3 chains. This is inferred from the coexistence of 
long and short range order, the field induced short range 
order, and the contrasting temperature behavior of the 
Q = and Q = X scattering. When the field pins the 
a spins into a long range ordered structure, the mani- 
fold of degenerate states is restricted considerably and so 
this also causes the (3 spins to become more correlated. 
The correlations in the (3 chains cause the scattering at 
the Q = X positions. Although the evolution of the or- 
dered Q = phase is not identical in all hysteresis loops 
the a chains always reach the same fully polarized state, 
independent of the state of the (3 chains. 

In both Dy2Ti2C>7 and Ho2Ti207 the scattering at the 
Q = X positions appears immediately, on application of 
a field, as a localized feature which then sharpens. This 
means that the interchain correlations begin to form as 
soon as the chains appear. This observation may be com- 
pared to those of Hiroi et al. who made a calorimetric 
study of Dy 2 Ti 2 7 in a field applied along [110}iS. On 
the basis of these results they argued that the chains 
behaved effectively as an assembly of one dimensional fer- 
romagnetic Ising chains embedded in a three dimensional 
crystal. We have directly measured the intrachain corre- 
lation length: at 1.5 T and 0.05 K it is of order 800 A for 
Dy2Ti2 07. The very long but finite intrachain correla- 
tion length is indeed a characteristic property of the one 
dimensional Ising ferromagnet. Our results clearly show 
that there are also correlations between the chains. At 
1.5 T and 0.05 K the interchain correlation length is of 
order 15 A for Dy2Ti207. The actual value of the corre- 
lation length depends on the history of the sample. The 
quoted values were at 0.05 K, longer correlation lengths 
are observed in the hysteresis loops at 0.3 K. 

Hiroi et al. pointed out that these chains are arranged 
in a triangular array, which frustrates the coupling be- 
tween them, preventing the development of true long 
range order. The triangular arrangement is not perfectly 
equilateral, so that in theory a transition to Q = X or- 
der is observed in the (3 chains. Our results also show 
clearly that the spin system does not achieve the long 
range order predicted at this point in (B, T) space 24 ' 25 . 

Indirectly these experiments have probed the dynam- 
ical processes operating in the spin ice materials in this 



very low temperature regime. The unusual situation of 
having the f3 chains perpendicular to the field and yet 
coupled to it by their interaction with the a chains gives 
particular insight into the dynamical regimes. These 
measurements show that there are dynamical processes in 
these materials at this extremely low temperature which 
have not yet been explained or detected in microscopic 
detail. The dynamics are extremely slow as these exper- 
iments operate at the longest timescales and yet equilib- 
rium is not maintained. Because we are below T = J e g 
the density of single spin flips must be exponentially 
small so they must be correlated dynamical processes. 

At 0.05 K both materials appear to be frozen: the cor- 
relations in the f3 chains only evolve when dynamics are 
restored by the field, via the a chains. As the tempera- 
ture is raised dynamics are increased. The boundary of 
the fully frozen regime was not clearly established. In all 
the experiments the histories followed placed the magnet 
in a strong field at very low temperature, with relatively 
disordered j3 chains (i.e. there was considerable potential 
for further order to develop in the (3 chains, as occurred 
subsequently) . In both materials this lead to the observa- 
tion that the (3 chains can be freed from the a chains and 
exhibit their own independent behavior. This was most 
clearly shown in Ho2Ti2 0y where the rising temperature 
lead to a startling increase in sharpness and intensity of 
the Q = X features. Similar observations were made 
in Dy2Ti207 where the degree of correlation achieved in 
the (3 chains was greater in a hysteresis loop at 0.3 K 
than at 0.05 K. As far as they are known, the measured 
dynamics of spin ices are consistent with thisM 4 * 4 ^* 4 ^. 
Dynamical processes appear to persist to lower tempera- 
tures in Dy2Ti2 07 which is consistent with susceptibility 
measurements. 

During all the experiments the spin ices appear to seek 
a long range ordered groundstate, but never realize one. 
This is shown by the ever increasing degree of correla- 
tion in the (3 chains when dynamics are restored and 
the a chains are controlled by the field. The increas- 
ing correlations are characteristic of the formation of the 
Q = X structure. The transverse ordering must be con- 
trolled by long range interactions, as the chains are not 
coupled by nearest neighbor interactions. It remains un- 
clear whether the equilibrium state is the fully ordered 
Q = X structure predicted by Yoshida— rather than the 
extended-range correlated state actually observed. The 
work of Harris et al. also should be examined 1 .. They 
concluded that the Q — X structure was a high tem- 
perature modification of the Q = structure. Here it 
is shown that because of the non-equilibrium nature of 
the low temperature regime, the magnetic structure ob- 
served is very much controlled by the dynamics and his- 
tory of the system. The question of the existence and 
spontaneous selection of a groundstate in spin ice remains 
open. It is surprising that after these extensive travels 
in (B, T) space, the materials seem to be arrested on the 
brink of order. Even at T w 1 K where the susceptibility 
is large and there are considerable dynamics available in 
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Ho2Ti207 it is unable to locate an ordered groundstate. 
It would be interesting to establish how the material be- 
comes jammed in this state and why the dynamics are 
unable to extend the correlations to long range order. 

At low temperatures (T < 0.5 K) the magnetization 
of the a chains depends on the rate of sweeping the 
field. In all hysteresis loops observed in this regime (ex- 
cept one with a slow field sweep rate in Dy2Ti2 0y) the 
a chain magnetization develops discontinuously through 
a series of jumps. In Ho 2 Ti 2 07, where repeated cycles 
were made, the pristine order parameter curve lies out- 
side of the main hysteresis loop which shows that the zero 
field cooled spin ice state is harder to magnetize than 
the field-cycled state. Neutron scattering cannot iden- 
tify if the order parameter jumps involve symmetry sus- 
taining first order phase transitions of the type discussed 
in Ref. ITU or they are a purely kinetic effect, involving 
the abrupt growth of macroscopic domains. If the lat- 
ter, the behavior of spin ice in this regard is reminiscent 
of random anisotropy Ising magnets such as amorphous 
Dyo.41Cuo.5gSl where similar magnetization jumps (like 
a giant Barkhausen noise) are observed. 

Field-induced magnetization steps with sweep rate de- 
pendence have also been observed in the systems Gd5Ge4 
and Lao.6Cao.4Mno.96Gao.o40;j£i. Both those systems 
have a field induced crossover from an antifcrromagnetic 
to a ferromagnetic groundstate, accompanied by a struc- 
tural transition. The steps were attributed to the burst- 
like growth of ferromagnetic regions in an antiferromag- 
netic matrix. The sweep rate dependence is believed to 
be due to the difficulty of accommodating the structural 
intergrowths at faster field sweep rates. Such behavior 
can also be attributed to the existence of two competing 
groundstates^. Here we have no structural phase transi- 
tion but there are two distinctive states at either end of 
the field sweep, the entropically favored spin ice state in 
zero field and the energetically favored Q = X state in 
high field. The situation with the field applied along [001] 
is analogous but with a different field induced ground- 
state. We suggest that these transitions might be viewed 
in the same way. Long range order grows in bursts out of 
the disordered spin system. In the martensitic transitions 
the bursts are due to the fact that system is well below 
the energy scale of the electrostatic interactions which 
maintain the structure and so the dynamics of the struc- 
tural rearrangement are extremely slow. Here the system 
is well below the energy scale of spin-spin interactions and 
the establishment of the ice rules regime severely restricts 
dynamics by preventing simple spin flips. Consequently 
all the transitions are being controlled by the balance of 



spin-spin interactions and magnetic field, with the small 
effective temperature promoting metastability. 



VII. CONCLUSIONS 

We have investigated the field-induced order in the spin 
ices Ho2Ti 2 07 and Dy 2 Ti 2 07. The degeneracy removal 
scheme depends upon the direction along which the field 
is applied. As expected, the application of a field along 
[001] leads to the selection of a long range ordered state of 
the Q = type. Application of the field along [110] leads 
to the separation of the system into a and j3 chains. It has 
been shown that the (3 chains have a long but finite intra- 
chain correlation length and only a very short interchain 
correlation length, and that the degree of correlation is 
history-dependent. A pathway to spontaneous selection 
of a long range ordered groundstate was not found, de- 
spite partial constraint of the system by the field. 

A prevalent feature of the experiments was the strong 
history dependence, usually manifested as hysteresis on 
cycling the field. It is concluded that as the experiments 
operate on such a long timescale, the dynamics of a spin 
ice are extremely slow in this temperature range. These 
slow dynamics mean that metastable states and coexist- 
ing long and short range order could be observed in both 
orientations. It was also shown that there is a change of 
dynamical regime from one in which the spin ice correla- 
tions are established but dynamical processes operate, to 
one in which the spin ice is much more strongly frozen. 
For measurements such as these this occurs at T « 0.5 
K for Ho 2 Ti 2 07 and T w 0.3 K for Dy 2 Ti 2 7 . The de- 
termination of the microscopic details of these dynamics 
remains a difficult problem. 
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